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ABSTRACT 

Context. Thirteen years after the discovery of the first afterglows, the nature of dark gamma-ray bursts (GRB) still eludes explanation: 

while each long-duration GRB typically has an X-ray afterglow, optical/NIR emission is only seen for 40-60% of them. 

Aims. Here we use the afterglow detection statistics of the systematic follow-up observations performed with GROND since mid-2007 

in order to derive the fraction of "dark bursts" according to different methods, and to distinguish between various scenarios for "dark 

bursts". 

Methods. Observations were performed with the 7-channel "Gamma-Ray Optical and Near-infrared Detector" (GROND) at the 2.2m 
MPI/ESO telescope. We used the afterglow detection rate in dependence on the delay time between GRB and the first GROND expo- 
sure. 

Results. For long-duration Swift bursts with a detected X-ray afterglow, we achieve a 90% (35/39) detection rate of optical/NIR after- 
glows whenever our observations started within less than 240 min after the burst. Complementing our GROND data with Swift/XRT 
spectra we construct broad-band spectral energy distributions and derive rest-frame extinctions. 

Conclusions. We detect 25-40% "dark bursts", depending on the definition used. The faint optical afterglow emission of "dark bursts" 
is mainly due to a combination of two contributing factors: (i) moderate intrinsic extinction at moderate redshifts, and (ii) about 22% 
of "dark" bursts at redshift >5. 

Key words. Gamma rays: bursts - Techniques: photometric 



1. Introduction 

Long-duration gamma-ray bursts (GRBs) are the high-energy 
signatures of the death of some massive stars, and they emit 
the bulk of their radiation in the a;300-800keV range. For 
understanding the physics of the GRB explosion, the impact 
of GRBs on their surrounding, as well as the implications of 
GRBs on early star formation and cosmology, it is crucial to 
observe their afterglow emission. While X-ray afterglows are 
detected basically for each burst, the low detection rate of op- 
tical/NIR af terglows has been a long standing problem in the 



GRB field dPaczvfiskiL Il998t IGroot et al l 1 1 998b iKlose et al 



200(1 iDiorgovski etajTl200lb iFvnbo et al.Ll200lb ILazzati etal 



2002t iPiro et all 120021: IKlose et a l.. 2003; iL^mejUJi [2006; 
Jaunsen et al. . 2008t iTanvir et all l2008l). an d it was first dis- 
cussed systematically in IFvnbo et all (l200lT) and ILazzati et al.l 
(12001 . 

The Swift satellite dGehrels et all 120041) was designed to 
slew to GRB locations rapidly and provide positions with arcsec- 
accuracy through ob servations of the aft erglows with the X- 
ray telescope (X RT. iBurrows etail 120 05) and the UV-optical 
telescope (UVOT. fRoming et all 12005). While the detection of 
the X-ray afterglows has dramatically facilitated the discov- 
ery of many new phenomena and increased the efficiency of 



ground-based follow-up observations, the UVOT detection rate 
of afterglows is, original ly somewhat surprisingly, only -40% 
(iRoming & Masonl 120061) . 

The first burst in the afterglow era for which no optical (o r 
NIR) afterglow was found was GRB 970828 (IGroot et all [i998). 
Originally, those GRBs with X-ray b ut without optica l after- 
glows had been coined "dark GRBs" dFvnbo et all 1200 lb . Later, 
this nomenclature was made more specific by adding a time and 
brightness limit, e.g. fainter than R ~ 23 mag within 12 hrs of 
the burst. As the next step, the b asic prediction of the fireball 
scenario (M eszaros &ReesLll997l) was utilized and the optical- 
to-X-ray spectral index fiox (slope between the fluxes in the R- 
band and at 1 keV) wa s used to define dark bursts (fiox < 0.5; 
iJakobsson et all 120041) . or the X-ray spectral and decay indices 



were used alternatively to e xtrapolate the X-ray flux to the opti- 
cal regime (Rol et afl, l2005l) . Ultimately, by using the X-ray flux 
and sp ectral (J3x) information from Swift, Ivan der Horst et afl 



ana sp < 

(2009) propose to define dark bursts as those for which fiox is 
shallower than (fix - 0.5). 

The darkness in the optical can have several origins (e.g. 
Fvnb o et all HOOI). The afterglow could (i) have an intrinsi- 
cally low luminosity, e.g. from an optically bright vs. optically 
dark dichotomy; (ii) be strongly absorbed by intervening ma- 
terial, either very locally around the GRB or along the line- 
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of-sight through the host galaxy, or (iii) be at high redshift 
(z > 6) so that Lya blanketing and absorption by the intergalactic 
medium would prohibit d etection in the frequently used R band 
dLamb & Reichartil2000t) . 

Before the Swift era, an analysis of a sub-sample of GRBs, 
namely those with particularly accurate positions provided with 
the Soft X-ray Camera on HETE-2, reve aled optical afterglows 
for 10 out of 11 GRBs dVillasenoretalll2004h . This suggested 
that the rapid availability of precise positions plays a major 
role in the identification of afterglows, and that the majority 
of dark GRBs are neither at high redshift nor strongly ab- 
sorbed, but just dim, and/or rapidly decaying. However, Swift 
observations have provided 500 GRB localizations at the few 
arcsec level within minutes of the trigger, and the fraction of 
non-detected afterglows is still abo ut ~30-40%: UVO T de- 
tects about 40% of the afterglows dRoming et all [2009), and 
ground-based observations recover another 20-30% (see e.g. 
http :/ w w w. mpe . mpg . de/~j eg/ grbgen .html) . 

This would imply that the accuracy and timeliness of a GRB 
position is not the dominating factor. On the othe r hand, given 
the larger mean redshift of G RBs in the Swift era dBerger et all 
l2005t|jakobsson et alll2006h . one could also argue that the mean 
flux of the afterglows is lower, and thus the effect of the better 
and faster Swift localizations on the afterglow detection rate is 
over-compensated by fainter afterglows with respect to the above 
HETE-2 sub-sample. 

Here, we use the afterglow detection rate of GROND, a de- 
dicated ground-based GRB follow-up instrument, to derive new 
constraints on the fraction of dark bursts. We complement the 
GROND data with Swift/XRT spectra to construct broad-band 
spectral energy distributions and to derive the rest-frame extinc- 
tion. Throughout this paper, we use the definition F v oc v~P. 



2. GROND observation statistics 

GROND, a simultane ous 7-channel optical/near-infrared imager 
dGreiner et alll2008al) mounted at the 2.2 m MPI/ESO telescope 
at La Silla (Chile), started operation in May 2007. GROND has 
been built as a dedicated GRB follow-up instrument and has ob- 
served basically every GRB visible from La Silla (weather al- 
lowing) since April 2008. 

GROND o bservations of GRBs within the first day are fully 
automated (see lGreiner et al.l2008al for more details). The fastest 
reaction on a Swift alert so far was 140 s (Swift trigger 353627); 
more typical times are 200-300 s for night-time bursts. The two 
dominant contributors for this delay are the read-out time of the 
interrupted exposure (particularly the Wide-Field Imager, one of 
two other instruments at the telescope), and the dome rotation 
speed. The distribution of delay times of GRBs with GROND 
is given in Tab. Q] (last column). Note that due to La Silla be- 
ing located at similar geographic longitude as the South Atlantic 
Anomaly (where most gamma-ray detectors need to switch off 
their high-voltage), the fraction of night-time bursts is nearly a 
factor two lower than for observatories at other geographic lon- 
gitudes. 

It is obvious that the later after a GRB trigger the GROND 
observation starts, the lower is the fraction of afterglow detec- 
tions (Tab. Q]). This is expected and readily explained by the 
rapid fading of GRB afterglows and the limiting sensitivity of 
the instrument/telescope in use. What is surprising, however, is 
the high detection rate in the first two time bins - this will be 
discussed in the following. The relatively high detection rate at 
>24 hrs is biased by the Fermi/LAT bursts which are, on av- 



erage, more energetic with correspondingly brighter afterglows 
(McB reen etaUl201Ch . 



Table 1. GROND afterglow detection fraction of long-duration 
bursts as a function of time delay of the start of the observation 
after the GRB trigger. 



Delay (hrs) 


detected vs. observed (fraction) 


fraction of total 


<0.5 


20/22 (91%) 


17% 


0.5^1 


15/17 (88%) 


14% 


4-8 


10/21 (48%) 


16% 


8-16 


22 /36 (61%) 


28% 


16-24 


13/22 (59%) 


17% 


> 24 


5 /10 (50%) 


8% 



Notes. Based on a total of 128 bursts observed between 070802 and 
(inclusive) 100331. 



Table [2] contains all those bursts which have been observed 
with GROND within less than 240 min (4 hrs) after the Swift 
trigger and which have XRT-detected afterglows (until 3 1 March 
2010). First we note that (within this sample) there is a clear 
difference between the detection likelihood for long- and short- 
duration bursts. While we detect only 1 out of 4 short bursts, 
it is 35 out of 39 for long-duration bursts. This resembles the 
well-known fact that short-duration bursts have generally much 
fainter afterglows than long-duration bursts. Second, we note 
that a further down-selection of long-duration bursts at high 
Galactic latitude does not change our detection rate. This is due 
to GROND's capability of imaging in all bands simultaneously, 
in particular also in the near-infrared. Third, we emphasize that 
this sample lacks any selection effects beyond requiring an X-ray 
afterglow. 

3. Fitting broad-band spectral energy distributions 

3.1. Overall properties 

GROND and SwiftPiKT data have been reduced in the standard 
manne r using pyraf/IRAF dTodvL Il993t iKupcii Yoldas et ail 
I2008bl) for the optical/NIR data and the XRT pipeline for the 
X-ray data. The optical/NIR imaging was calibrated against the 
primary SDSS standard star network, or cataloged magnitudes 
of field stars from the SDSS in the case of g'r'i'z' observations 
or the 2MASS catalog for JHKs imaging. This results in typi- 
cal absolute accuracies of +0.03 mag in g'r'i'z' and ±0.05 mag 
in JHKs ■ X-ray data were cleaned for time intervals of flaring 
activity and the early steep decay. The X-ray spectra were flux 
normalized to the epoch corresponding to the GR OND observa- 
tions using the XRT light curves from lEvans et a l] (l2007l l2009). 
This common reference time was selected to be after the early 
optical rise observed in some light curves and is different for 
each burst, but typically during the first few hours after the Swift 
trigger. 

We then combined X RT and Galactic foreground extinc- 
tion dSchlegel et all fl998h corrected GROND data to establish 
a broad-band spectral energy distribution (SED) - see Fig.Q] 
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Table 2. GRBs observed with GROND within 240 min after the Swift trigger (including 39 long- and 4 short-duration bursts). 
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Notes. Columns 2-10 are the duration (D) classification of the GRB according to the canonical long (L) / short (S) scheme 1 Kouveliotou et all 
Il993h . Galactic Latitude B, the delay between the GRB trigger and the start of the GROND observation (which is the sum of the delay between 
the burst and the arrival of the Swift notification at the GROND computer, and that until the GROND start), the duration (Dur) of the GROND 
observation during the first night, the afterglow (AG) detection by GROND (G) or others (O), references for the afterglow observations, the 
redshift z with reference, and special comments (UVOT detection or X-ray brightness). The last four columns are the best fit spectral slopes in the 
optical/NIR (Bo) and X-ray (J3 X ) band, the rest-frame extinction A v and absorption Nh from the combined GROND/Svw/f-XRT spectral fit. Bo was 
generally fixed to B x (-0.5) in the fit, and thus has the same error as B x - 

" References for prev iously reported GROND a fterglow detections as w ell as those from othe r gr oups for those GRBs whic h have not been de- 
tected b y G ROND: (l)lKruhler et alj <2008d). (2)ICovino et al l (120081), (3)lKriihler et al.1 d2009 st), (4)lGreiner et al.1 d2009bt). (5)lKupcu Yoldas et al.l 
( 2008 ^). (6) | Rossi et al.l d2008 j) . (7)lciemens et alJ i2008ah.lGuidorzi et al.1 d2009h . ( 8)lKruhler et alj i2008ah . (9) lRvkoff & Ruiopakard J2 008). 
(10) iKriihler et alj j2008bh , (1 1) iFilgas et alj j2008ah . (12) iRossi et al l d2008bl) . (1 3) IClemens et alj j2008bh. (14) IClemens et all ( 12008^)7(15) 
Kruhl er et al.l d2009tJ). (16)lKriihler et all i2008eh. (17)lKruhler et all i2008d|). (18)lRossi et al.l d20083)nG reiner et all d2008bl). (19))Rossi etaD 
d2008dh, (20)lDella Valle et al.1 d2008l).'(2 1 )lYuan et all d2010D, (22)IClemens et"5] d2008dh. (23)Il6w et al.1 d2008h, (24)|Updike et alj d2008h. (25) 
I Afonso et al J d2008ah. (26) I Afonso etaD d2008bh. (27)lAfonso et al.l d2009h: iGendre et alJ d2010h. (28) ICenko et al.l d2009 j). (29) lUpdike et ; ~ 
d2009bh. (30)IOlivares et al.1 d2009ar).731)lRossi etall d2009aT ) , (32)IUpdike et alj d2009d)~(33)IUpdike et al.1 d20093K ( 34)IOlivares et al.1 (120091: 



d2009eh . (39) IFilgas et al.l d2009cl) . (40) 



(35) Malesani et al. (2009), (36) Filsas et al. (2009a), (37) Filsas et al. ( 2009b), (38) Updike et al. 
iLevan et al.1 ( 120101) . (41) lNicuesa Guelbenzu et al.1 d2010l) : lTanvir et alj d2010l) , (42) lAfonso et alj d2010l) 

2) A "p" after the red shift value indicates a photo metr ic redshift es t imate . Reference s for t he r edshifts are: (1) iProchaska et all d2007al): 
Eliasdottie retaLl J2009]), ( 2) iProchaska et al.1 d2007bl), ( 3)lLedoux et all d2007[) ; lFox et all ^003), (4)lGreiner et all d2009bh. (5)Uakobsson et atl 
(2008 j), (6) iFvnbo et all d2009bl). (7) iMalesani et all d2008l): iGuidorzi et alJ O009I) . (8) PrhQne et all d2008al). (9 ) lTh6ne et all d2008b[). (10 ) 
IVreeswiik et al.1 d2008l), (11) IFilgas et all d2008bT)~ (12) Uakobsson et all d2008bl). (13) Ijakobsson et al l d2008cl) . (14) IFvnbo et all |2008 
(15) iPerlev et all d2008t>. ( 16) iThone et al] d2008ch (17) Ijakobsson et al.1 d2008d). ( 1 8) iGreiner et all d2009£<) . (19 ) iBerger et alJ d2008 
(20) iD'Avanzo et al.1 d200g) (21) iD'Elia et al.1 d20ol). (22 ) iBerger & Rauchl d2008l). (2 3 ) ICucchiara et al.1 d2008h (24) lAfonso et al.1 d2008d 
Kruh ler et all d2qi0bl). (25) Ide Ugarte Postigo et all d2009al). (26) IChornock et all d2009h : Ide Ugarte Postigo "et all |2010l) . (27) ICucchiara et al! 
d2010bh:lTanvirl d201( " 
(32) IChen et"aiTd2009l 

3) The values for host A v and N H are given for redshift zero. 



(28) Thone et al. (2009), (29) de Ugarte Postigo et al. (2009b), (30) Jakobsson et al. (2009), (31) Fynbo et al. (2009a), 
, (33) IChornock et all d2009bl) . (34) ICucchiara et al.1 d2009bl) . (35) rFong et all d20irE 736) Vergani et alj d2010h . 
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These SEDs have been fit with two alternative models: (i) 
a single power law with free slope and normalization, plus free 
source-intrinsic extinction of SMC/LMC or MW-type (for the 
GROND data) and Galactic plus rest-frame equivalent neutral 
hydrogen column density (for the X-ray data) assuming solar 
abundance; or (ii) a broken power law where the break energy is 
left free but the difference in the two slopes is fixed to 0.5, and all 
other parameters are left free as above. Note that fixing the break 
does not imply that the optical/NIR SED is not fit; instead, it is 
fit together with the X-ray SED (the X-ray slope in general does 
not dominate the SED fit). These models provide a very good fit 
to all but one afterglow SED (Fig.Q]): For this exceptional case, 
GRB 080413B, the combined GROND/XRT SED requires a dif- 
ference >0.5 between the l ow and high energ y slopes of the syn- 
chrotron emission model (Filgas et al.L [20Tof> . Except for seven 
GRBs (080605, 080710, 080913, 081029, 081228, 090926B, 
091221) all GRBs are better fit with a break between the X- 
ray and optical/NIR. The best-fit parameters are listed in the 
last four columns of Tab. [2] With the exception of five bursts 
(GRBs 070802, 080210, 080605, 080805, 090102), the SEDs of 
the afterglows are consistent with being reddened with an SMC 
extinction law. 

3.2. Notes on individual GRBs 

070802 has a fiox = 0.5, but since the best-fit location of 
the spectral break is just bluewards of the g'-band, it still re- 
quires a substantial Ay ~ 1.2 mag. This burst is also an exam- 
ple of a bor der-line case: its figxjust misses the "dark burst" 
criterion of Ijakobsson et alJ (120041) . but due to its fix being 
sl ightly steeper than 1 .0, it qualifies as "dark burst" according 
to (van der Horst et all d2009l) . 

080218B: If we assume a break of 0.5 between the X-ray and 
optical/NIR band, as for the majority of the bursts in our sample, 
the GROND upper limits allow the following (z, Ay) pairs of e.g. 
(3.5,1.5) or (5,1.0) or (7,0.7) or (10,0.5) as explanations of the 
non-detection. 

080516: With the improved calibration, our best-fit redshift is 
z=3.6+0.6. Since GROND is not sensitive to redshifts smaller 
than »3, the above redshift is consistent with zero at the 90% 
confidence level. However, we use this best-fit redshift to derive 
Ay which therefore could be a lower limit. 
080805: In this case there are no simultaneous GROND and 
XRT observations of the afterglow. The flux normalization of the 
XRT spectrum was obtained by back-extrapolating the XRT af- 
terglow light-curve to the earlier time of GROND observatio ns. 
080915 has a very faint X-ray afterglow dOates et aUl2008l) . If 
we assume a break of 0.5 between the X-ray and optical/NIR 
band, as for the majority of the bursts in our sample, we predict 
r 1 ~ 23.2 mag with Ay = 0. This is brighter than our GROND de- 
tection limit (r' > 25.2 mag) for this burst. Thus, various (z, Ay) 
pairs of e.g. (3.5,>0.5) or (5,>0.25) or (7,>0. 1) provide good ex- 
planations of the non-detection (Fig. [3]). 

090429B: Due to bad sky conditions, this burst was not seen in 
the first night, but detected in the second night at the 2cr level 
in the H band, but not in J. Since GROND is more sensitive 
in J as compared to H, this is consistent with a break between 
J and H, and thus with the photometric redshift determined by 
ICucchiara et"ail d2010bl) . 

090904B: The object is detected in all GROND filters except 
the g'-band which implies a secure upper limit on the redshift of 
z < 5. However, the g'-band non-detection is also well explained 
by th e large Galac t ic fo reground reddening (E(B - V) = 1.76 
mag lSchlegel etalj d!998h ). With no accurate redshift available, 



the best-fit extinction of Ay =2 mag at z — is an upper limit. 
However, the large foreground extinction of Ay ~ 5 mag can be 
expected to come with a large (systematic) error, which propa- 
gates also to our fit values; thus we cannot distinguish between 
Galactic foreground and host extinction. 

090926B: While GROND observations of this GRB started 
quickly, clouds prevented continued observations after a few 
minutes. Further observations were done starting 6.2 hrs af- 
ter the burst, and additionally on 14 Feb 2010. The latter ob- 
servation confi r med th e fading of the candidate counterpart of 
Male sani et alJ d2009l) . and the association to the host galaxy 
as well as the red shift determination via host spectroscopy 
dFvnbo et al.Ll2009al) . 

091221: No spectroscopic redshift is available. The afterglow 
is detected in all GROND filters, and the 3cr upper limit on 
the redshift is z<3.3. For a putative redshift of 2, we obtain 
A v = 0.21+° |3 mag and JV H = 1!| 4 10 21 cirT 2 . 
1001 1 7A: Th e redsh ift of z=0.92 is containe d in a summary ta- 
ble in lBergerl d2010) and referenced there as iFong et alJ fcOlO. 
in prep.). 

100205A: Due to cloud coverage, GROND observations started 
only 2.3 hr after the burst, and then were severely affected by 
passing cirrus. Correspondingly, no detection was achieved, and 
the upper limits obtained by GROND were about 2 magni- 
tudes worse than the nominal s ensitiv ity under normal condi- 
tions (Nic uesa Guelbenzu et all feOlOh . At about the same time 
(2.6 hr after the G RB), Gemini observations revealed a near- 
infrared counterpart iTanvir et"ai] d2010l) with a very red colou r 
of H-K= 1.6 + 0.5 mag (AB system - ) dCucchiara et all l2010a). 
If this colour is due to Lyman-Q' absorption within the H fil- 
ter bandpass, then this would imply a redshift of 11 < z ~ 
13.5. Lower-redshift solutions with significant local (host) ex- 
tinction, unusual afterglow colours, or substantial contributions 
from an underlying host galaxy have not been excluded though 
dCucchiara et alll2010al) . 



4. Discussion 

4.1. The nature of dark bursts 

In order to determine which of our burs ts would have be e n clas- 
sified as dark bursts according to the Jakobsson et al. (2004) 
definition, we also estimate /fox (Fig. |2]l using the measured 
7?-band brightness/limit and the X-ray flux at 3 keV (both ob- 
server frame). We find a fraction of 26% (the nine GRBs 
080218B, 080516, 080805, 080913, 080915, 090102, 090429B, 
090812, 090926B) that would be considered "dark". Another 
five bursts (070802, 080129, 080413B, 080605, 090814) have, 
within errors, all /3px = -5. According to the classification of 
Ivan der Horst et all (120091) . the following nine bursts are "dark": 
070802, 080210, 080218B, 080516, 080805, 080913, 080915, 
090429B and 090904B. These fractions of 25-40% are fully 
consistent with the hitherto known fraction of dark bursts of 
25%^12% (e.g. lFvnbo etafll2009bl) . 

There is one other problem with all of the "dark burst" 
definitions, namely, the temporal evolution. Usually, a certain 
time interval is selected, based on the availability of optical and 
X-ray measurements, and one of the "darkness" tests applied. 
However, given the increasing sample of bursts with not just 
good X-ray but also good, simultaneous optical coverage, the 
complexity of light curves as compared to the basic fireball sce- 
nario becomes more and more evident. We have seen all kinds of 
different light curve behaviour, with either the optical or the X- 
rays decaying faster; sometimes this even changes in one burst. 
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Fig. 1. Spectral energy distributions of 
a selection of the long-duration GRBs 
with redshift in our sample, labeled on 
the right side. The X-ray spectrum has 
been extracted from early times, but af- 
ter the rapid fading and excluding flares, 
and then shifted in mid-time to the 
GROND mid-time (see text for details). 
The dashed lines indicate the best-fit 
model for each burst, with the break be- 
tween the X-ray and optic al/NIR always 
treated as a free parameter. Curvature in 
the X-rays is due to Galactic plus host- 
intrinsic absorption, and curvature in the 
observed optical/NIR range due to host- 
intrinsic extinction (data have been cor- 
rected for Galactic foreground extinc- 
tion before fitting). For better visibility, 
the intensity scaling has been choosen to 
minimize overlap of the different SEDs, 
and thus is completely arbitrary. 



As a consequence, there is nothing like a single /fox per burst, 
but instead each burst shows a more or less pronounced evolu- 
tion of /fox- As an example, the behaviour of GRB 080413B is 
shown in Fig. [2] 

The main results of the SED fitting with respect to the 
"dark" burst isssue can be summarized as follows (the follow- 



ing properties are n on-exc lusive, and we use the definition of 
Ivan der Horst et al] d2009f) in the following): (1) Four out of 
the nine dark bursts are faint due to non-zero, but moderate 
(A v ~ 0.2 - 1.5), extinction: 070802, 080210, 080516, 080805. 
The measured extinction in many cases appears enhanced to the 
observer due to a moderate redshift of the burst (see Fig. [3]). (2) 



5 



J. Greiner et al.: The nature of "dark" gamma-ray bursts 




10 1 10* 10 s 10* 



K,x My] 



Fig. 2. The distribution of our sample GRBs in th e plane of 
fl-band flux over X-ray flux iJakobsson et al.1 (120041) . For GRB 
0804 13B we show as a thick solid line the temporal evolution of 
ySox with time during the first 2 days of the afterglow evolution 
(from top to bottom left). 



One burst (090904B) is behind an additional Ay ~ 2 mag dust on 
top of the nominal Ay ~ 5 Galactic foreground which could ei- 
ther be due to patchiness of the foreground or due to host galaxy 
extinction. (3) Two bursts are faint (in r') due to high redshift 
(z > 5), namely GRB 080913 and 090429B. This corresponds to 
a fraction of 22% +8% of the dark bursts. (4) The remaining two 
of our "dark" GRBs, 080218B and 080915, are either at large 
redshift or moderate Ay, or both (see notes above) - so they be- 
long to one of the above two groups (1) or (3). (5) Even bursts 
with good evidence for a spectral break and fiox = 0.5 do re- 
quire extinction in some cases, in particular if the break is near 
the optical (rather than X-ray) range. 

4.2. Extinction 

The rest-frame Ay could be derived for 33 bursts of our sam- 
ple (TabJSJljandjts distribution is compared against that of 
the iKann et ail (1201 Oh collection (Fig. @}, made up of all op- 
tically bright GRBs with photometry available from the litera- 
ture. We find substantially larger extinction values, in particular 
about twice as many bursts (~25% vs. 12%) with Ay ~ 0.5 mag, 
and for the first time a significant fraction (10%) of bursts with 
Ay > 1 mag from a direct optical/NIR afterglow SED redden- 
ing measurement. The difference is easily understo od as a se- 
lection effect (predominance of bright bursts) of the IKann et aD 
d2010l) sample. Note, however, that this result is not due to the 
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Fig. 3. The effect of the combination of various values of in- 
trinsic extinction Ay (lines with labels) and redshift, produc- 
ing an effective dimming in the R band as given on the y-axis. 
The solid lines have been computed assuming a dust extinc- 
tion curve as described in iReichartl (l999) but with no 2175 A 
bump and an opacity due to intergalactic hydrogen approximated 
by t(HI) = 2.6 x (1 + z)" for 700 nm/(l + z) < 121.6 nm 
dValageas. Schaeffer & SilMll999h . Dots represent the GRBs of 
our sample for which we have an Ay measurement from the 
GROND SED. GRBs with effective R-band reduction of >0.5 
mag are labeled; GRBs 080218B and 080915 lie above the cor- 
respondingly labeled dashed lines which correspond to the pairs 
of (z, Ay) values described in the text. Moderate Ay at moder- 
ate redshifts can easily produce a dimming of 1-3 mag in the R 
band. 



way in which the extinction models are fit to the data - our re- 
sult is derived with the canonical SMC, LMC or MW extin ction 
laws, in the same way as was done in IKann et alJ (l2010h . We 
note in passing that using a different fitting approach (a 'Drude' 
model) was reporte d to yield a factor 2-5 larger visual extinction 
(iLiangetallfioiOh . 

It is also interesting to note that the bursts with the largest 
measured Ay are those which tend to prefer extinction laws that 
include a 2175A feature. Unfortunately, statistics is poor, so one 
has to await for future, larger samples to verify whether or not 
this is a generic trend. 

Despite this larger fraction of moderate to large extinction, 
there is still a 50% fraction of bursts with Ay = 0, with er- 
rors of only +0.02-0.05 mag (Tab.|2]i. This has also been some- 
what surprising for two reasons. First, in the canonical picture 
of a massive progenitor with a strong wind prior to explosion, at 
least some extinction from the eje cted material w ithin the wind 
of the progenitor is expected (e.g. Waters, 2010, and references 
therein). Second, even if the local environment of the GRB is 
relatively dust free, the GRB radiation has to pass through the 
host galaxy to reach us. Already the latter effect alone should 
create a mean ex tinction of order Ay = 0.05 - 0.1 mag for about 
30% of sources 

(UpdikeetaU,|2009a), 

and very few with high 
Ay; see Fig. [4] We now observe both, a larger fraction of Ay = 
objects as well as a larger fraction of Ay > 1 objects as com- 
pared to this, admittedly simplistic, expectation. The reason for 
this discrepancy is presently unknown, but could be related to 
the possibility that the dust geometry is not homogeneous, but 
clumped. The relative overabundance of Ay = and Ay > 1 
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Fig. 4. The distribution of mea sured Ay from ou r sample with 
redshift as compared to that of iKann et al.l d2010l) . The dashed 
line is the theoretical distribution, normalized to the same num- 
ber of objects in the sample, of a Monte Carlo sim ulation of ran- 
dom s ightlines through an evolving galaxy model (lUpdike et al.L 
2009a). A KS-test returns a rejection probability of 54% for the 
null hypothesis that the two distributions are drawn from the 
same sample. 



objects would then be determined by the covering factor and the 
dust column through the clumps. 

It is inte resting to note that the Ay distribution of type II-P 
supernovae dSmartt et al.Ll2009h (which comprises about 60% of 
all core-collapse SN), an d that of stripped-envelope supernovae 
dRichardson et all 12006) (which comprises 3 lib, 11 lb and 13 
Ic) derived from a comparably sized sample as our GRBs, are 
consistent with being drawn from the same underlying distribu- 
tion as our Ay distribution of GRBs (though many of their Ay 
values have errors consistent with Ay = at the lcr level). This 
suggests that we sample similar environments for both popula- 
tions, hinting at similar progenitor systems: all these supernovae 
tend to be associated with star-forming regions, so they can also 
be expected to be significantly extinguished in their host galax- 
ies. There is similar evidence that the distribution of Wolf-Rayet 
stars is consistent with the theoretical picture that type-Ic SN re- 
sult from prog enitors that have been stripped of a larger part of 
their envelope dLeloudas et all 1201 Oh . 

4.3. Gas-to-dust ratio 

Effective hydrogen absorption in excess of the Galactic fore- 
ground absorption has taken a long time to be detected signifi- 
cantly in GRB afterglow spectra. Origi nally not detected at all in 
the full sample of BeppoS AX bursts (De Pasauale et al.Ll2003l) . a 
re-analysis of the brightest 13 X-ray afterglows revealed statisti- 
cally s ignificant absorptio n in excess of the Galactic one for two 
bursts (IStrattaetal.Ll2004 . Already 8 bursts of 17 observed with 
Chan dra or XMM - Newto n up to Oct. 2004 show excess absorp- 
tion dGendre et al.L I2006I) . In a systematic study of 93 promptly 
observed Swift GRBs with known redshift (up to May 2009), 85 
show evidence of intrinsic X-ray absorption at the host galaxy 
site dCampana et all 120101) . Similarly, in our sample, we detect 
excess absorption in 26 out of 33 cases. This difference in the ex- 
cess absorption detection rate is primarily related to the quality 



Fig. 5. Relation, or lack thereof, between best-fit rest-frame 
extinction and rest-frame neutral hydrogen absorption of 
our sample. The lines are the relation for our Galaxy 
dPredehl & Schmitt L fl995l) . as well as lOx and lOOx larger N H 
(from bottom to top). GRBs 080129 and 080516 are missing as 
they fall outside the plot with their large A^h values. 



of the X-ray spectra, which pie-Swift was typically taken 8-12 
hours post trigger, by which time the signal-to-noise ratio was 
insufficiently high to accurately measure any intrinsic absorp- 
tion. 

Figure [5] shows a comparison of the effective neutral hydro- 
gen absorption A^h and visual extinction Ay, both in the rest 
frame of the corr esponding burst. A s has been noted frequently 
in the past (e.g. iGalama & Wiiersl 1200 U IStratta et all 120041: 
Scha dv et all 120101) . the A^-to-Ay ratio is far from being sim- 
ilar among different bursts, and substantially larger than in our 
Galaxy. Note, however, that here, as has been usual in previous 
cases, solar metallicity has been assumed in deriving Nh- Since 
the observed curvature in the X-ray spectra is predominantly, 
but not exclusively, due to absorption by oxygen, the derived ef- 
fective A^h is inversely proportional to the metallicity (or better 
O/H ratio) of the burst environment. Since this metallicity has 
been observed ( in other GRBs) to be about 1/5 solar (though ex- 
tremes of solar (Procha ska et all l2009t)Eliasd6ttier et all 12009) 
up to super- s olar (Savaglio et al., 2010) and nearly 1/100 solar 
dD'Elia et all 120071: iRau et all 1201 Oh do also occur), the effec- 
tive Nh would likely be even larger than shown in Fig. [5] if 
the proper line-of-sight metallicity were to be used (if it were 
known). 

In contrast to our Ay distribution, the distribution of A^h.x 
from the complete S wift sample lacks a su bstantial fraction of 
zero colum n density (jCa mpana et al., 2010). This has been ex- 
plained by ICampana et alj (|2010) as evidence that the bursts 
originate within high-density regions of their hosts, since a ran- 
dom distribution in a galaxy like ours would predict a siz- 
able fraction (~30%) with no intrinsic absorption. By com- 
bining their sample with the Lyman-ff absorbers at z > 2 of 
Fvnb o et all d20 09b) they also find, similar to earlier reports (e.g. 
Wats on et all [2007), that the bulk of GRBs have column densi- 
ties in X-rays which are a factor ~10 higher than in the optical 
(A^hi), which they explain by ionization of hydrogen by the high 
energy flux of the GRB. Since this ratio is roughly similar to 
that of A^h.x vs. Ay (Fig. [5]), one could think that Ay and the 
Lyman-a absorption are correlated. However, the three bursts 
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with reported N m in lFvnbo etal] d2009bl) (070802, 071031 and 
080804) do not show any correl ation, similar to the 6 bursts from 
the UVOT sample published bv lSchadv etail d2010h . 

4.4. Redshift distribution 

The redshift distribution of our sample is 92% complete (also 
one of the 4 GRBs not detected by GROND has a redshift). We 
re-iterate that the only selection criterion was the detection of an 
X-ray afterglow, and do not see a bias introduced by the require- 
ment of a rapid GROND observation (or equivalently an occur- 
rence during Chilean night time). A comparison to the distribu- 
tion of all known long-duration bursts (about 50% complete; Fig. 
[6]l reveals the former to have a flatter distribution, with a some- 
what higher number of z > 4 bursts. However, a KS-test shows 
that this is not statistically significant, and both distributions are 
consistent with being drawn from the same sample within lcr. 

The presence of GRBs 080913 and 090429B in our sample 
corresponds to a fraction of 5.5+2.8% of bursts at redshifts z > 
5, and the strict upper limit would be 12.8% if all three GRBs 
without redshift would be at z > 5. A larger sample size with a 
similarly good completeness level would be required to derive 
fractions with errors less than the present ~50% level. 
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This finding is in line with previous investigations 
elandri et all 12001 ICerA oltaLl l2009bt iPerlev et all 120091: 
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Fig. 6. Redshift distribution of our GROND sub-s ample (blue) 
vs. tha t of the hitherto most complete sample of (Fvnb o et al 
2009b) (bursts with upper limits have been omitted). 



dMelandri et all 120081: I Cenko et al. 

iFvnbo et all l2009bl) . In particular. [Cenko et all (l2009bl) used a 
similar approach as ours, namely a sample of 29 bursts for which 
follow-up observations with the robotic Palomar 60 inch tele- 
scope began within 1 hr after the burst trigger. They recovered 
80% of the optical afterglows, as compared to our 90% with a 
~60% fraction of redshifts (as compared to our 92%). A search 
for host galaxies was then performed for bursts without after- 
glow and/or redshift of this sample to as sess the degenerac y be- 
tween high-z/low-Ay vs. low-z/high-Ay dPerlev et all 120091) . In 
combination, and assuming further that the extinction measured 
along the line-of-sight to GRBs is proportional to the extinction 
of the host emission, these authors reach the similar conclusion 
that extinction is responsible in large part for the "dark bursts". 
Some of their Ay values, however, have been derived by assum- 
ing a SED slope in the optical, rather than measuring it. Given 
the prevalence of breaks (see above), this is risky. 

In contrast, we emphasize here, that we (1) measure the op- 
tical/NIR SED, and (2) make no assumptions on the host prop- 
erties. Thus, our sample of bursts is the first with properly mea- 
sured Ay values which neither requires a relative shifting of dif- 
ferent filter measurements, as often done in the past nor suffers 
from the small wavelength coverage. Our only assumption is 
that whenever the data prefer a break in the SED between X- 
rays and optical/NIR, we only allow a slope difference of 0.5. 
GROND observations start removing the bias of previous stud- 
ies to bright bursts, i.e. low-Ay and low-z, and start to detect 
and reliably measure these higher-Ay values directly. From in- 
dividual bursts not in this sample there is evid ence that bursts 
with e ven A v ~ 3 - 5 exist (e.g. GRB 080 607, iProchaska et all 
(2009); GRB 070306. Haunsen et al l (120081) ') - though these are 
admittedly rare occasions which require larger samples to in- 
clude them. 
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5. Conclusions 

The majority of afterglow SEDs show a spectral break between 
the X-ray and optical/NIR range that can be described with a 
slope difference of 0.5, consistent with the basic fireball sce- 
nario. This spectral break implies a /?-band flux of about 3-4 
mag fainter than obtained from an extrapolation of the X-ray 
spectrum. This effe ct is dealt with in all the recent definitions 
of "darkness" (e.g. Ped ersen et all 120061) . and our finding of a 
dominance of SED breaks is consistent with the large number of 
bursts seen at f3o < 1 . 

The faint optical afte rglow emission of "dark bu rsts", where 
we used the definition of Ivan der Horst et al. (2009), is due to a 
mixture of moderate intrinsic extinction at moderate redshifts, 
and a fraction of bursts at redshift >5 (about 25% of the dark 
bursts). 



Facilities: Max Planck:2.2m (GROND), Swift 
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